Exactness of the replica method in perturbation 



Hisamitsu MukaidaE 
Department of Physics, Saitama Medical University, 38 Mora hongo, 
Moroyama-cho, Iruma-gun, Saitama, 350-0495, Japan 

Yoshinori Sakamotclll 

Laboratory of Physics, College of Science and Technology, Nihon University, 
7-24-1, Narashino-dai, Funabashi-city, Chiba, 274-8501, Japan 
(Dated: 090112 vl) 

0\ ' 

' The replica method for a quenched disordered system is considered in a perturbative field theory. 

, Since correction in a finite-order perturbation is given in a polynomial of the replica number n, the 

■ zero-replica limit n — > is regarded as extracting the constant term from the polynomial, which 

, ' mathematically makes sense. The meaning of the extraction is clarified comparing with a direct 

S ' calculation. 



C/3 



> 



PACS numbers: 05.10.Cc, 64.60.Ae, 64.60.De, 64.70. Q- 



I. INTRODUCTION 



The replica method is widely used for studying quenched disordered systems It was developed for getting 
around difficulty involved in taking average E (•) over quenched disorder. For example, in the case of a free energy, 
we compute E Z" instead of E log Z, where Z is the partition function of the system computed under fixed quenched 
disorder, and use the formula 



'rt i ElogZ = lim i(EZ" - 1). 

o 
o 



Here, is realized by introducing n "replicants" identical with the original system, so that n is a positive integer, 
and thus the limit 7i ^ is ill-defined. Although there are specific models in which the replica method is proven to 
be mathematically rigorous [1, |3, 01, justification for the zero- replica limit is left beyond the scope of investigation in 
most literatures. Nevertheless, since this method gives rehable results to various quenched disordered systems, it has 
survived for more than three decades. In this paper, we wish to point out that the replica method can be justified 
within the framework of perturbation. 
, The main idea for the justification is that a physical quantity in the replicated system calculated in a finite-order 

■ perturbation is a (finite-degree) polynomial in n. The zero-replica limit can be regarded as the extraction of the 
Q,^ ' constant term from the polynomial. A similar idea is used by Brunet and Derrida in Ref.f^. 

■ Generally, if a function f{n) defined on positive integers is a polynomial with the degree m, the constant term is 
obtained solving a linear equation for all the coefficients of f{n) generated by /(I), ...,f{m + 1). We thus find that 
the constant term is extracted by 



m+l 



' n—l ^ ^ 

P.: 

It should be noted that Vo in ([T]) becomes ill-defined when m = oo, which implies that the constant term of a 
power series in n becomes ambiguous. For instance, consider the series J2T=o(^^)'^ i^^)'^'' / i"^^ ^ ^V- ' which apparently 
indicates that its constant term is one. On the other hand, the summation results in sin n7r/(n7r), which vanishes for 
all n = 1, 2, .... Such ambiguity causes probability densities where the replica method does not work 

In this paper, we show that the constant term in n of a physical quantity calculated in the replicated system equals 
the disordered average of the corresponding physical quantity in the original system within finite-order perturbation. 
It is also found that the n-dependent terms in the replicated system originate from the disorder correlation between 
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the corresponding physical quantity and the free energy in the original system, which are nothing to do with the 
disorder average we want. The limit n —^ can be interpreted as realization of extracting the constant term. 

Our argument is also applicable to perturbative renormalization group, which includes functional renormalization 
group in perturbation j^, Q . We will see that a beta function calculated in perturbation becomes a polynomial in n in 
the replicated system. The constant term in it is precisely equal to the beta function for the correlators characterizing 
probability density of quenched disorder. 

It is pointed out in the literature that a replicated field theory describing the random-field Ising model (RFIM) 
yields one-loop correction to coupling constants singular in the limit n — > 0. A novel replica limit is proposed in order 
to remove the singularity. Since it provides beta function allowing breakdown of dimensional reduction for the RFIM 
near the upper critical dimensions contrary to the conventional limit it is important to know which the right 
way is. Although the two replica limits can be compared with a direct computation not employing replicants [T^ . the 
ambiguity could not be resolved within the replica method so far because of missing mathematical basis of the limit. 
This issue will be treated as an example. 



II. EXPECTATION VALUE 



Let us consider a field theory consisting of a single-component field (j> and a quenched disorder v. We write the 
thermal expectation value of a physical quantity A[(j)\, a functional of 0, as 

where H[(f); v] is the Hamiltonian (times l/kT) of this model. We can put -ff [0; v] = for an arbitrary v without loss 
of generality since the quotient ^ is independent of H[Q;v\. 

The probability density for v is usually characterized by a given set of correlators among itself, which is denoted by 
u. The average and the cumulants for the quenched disorder are respectively described as E and k. E.g., k {v) = Ew, 
K (wi, V2) = E (vi V2) — EtiiEt;2- 

Our interest is to compute expectation values such as E (^[0])// [0.^,] by means of the replica method. Introducing 
n identical systems, we define the replica partition function 

/n 
]J l?0c,e-^°=i^['^°'^l. (3) 
a=l 

The replica Hamiltonian is defined as 

-'H[ci);u] = logEe"^"=i'f^['^°;"l, (4) 

where we have used the notation = [(pi, </)„). We restrict ourselves to the case where the cumulant expansion for 
the right-hand side terminates at some finite number: 



^ai — 1 



which excludes introducing infinite number of replica indices. The physical quantity in the replicated system corre- 
sponding to E is (^[0i])-H[0.u] ' where the subscript 1 is the replica index. It can be written as 



Now we evaluate the right-hand side in finite-order perturbation. Letting W — log Z, the cumulant expansion has 
the form of 



^==0 ■ V I J 



(7) 
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We write the unperturbed Hamiltonian as Hq [(j>; u] . Here, the expression Hq [(p; u] instead of Hq [(f); v] imphes that 
the disorder average for Hq is aheady taken. Perturbation is defined as ^ = [(/);«] — Hq[(I);u]. The entries in the 
cumulant ([7]) have the following perturbative series: 



M^ = logZo+ J2 



i-iy 



A[4>];V;..-V 



™ Ho\(i:u\ 



v-.-v 



(8) 



which are obviously independent of n. Here the semicolons in the angle brackets mean to take the connected part, 
e.g., {A;B) = (AB) — (A) (B), and the partition function for Hq is denoted by Zq. The summations in ([7]) and ([5]) 
are truncated at a finite number under a finite-order perturbation, so that is expressed in a polynomial 

in n. Let [Y]m be the perturbative expansion of an arbitrary quantity Y up to and including the m-th order in V. 
Since the constant term of ([7]) is given by the term with k = 0, we have 



(^[<^])h[0;.] • 



(9) 



It shows that the perturbative expansion for E 



is given by the constant term in ri of the perturbative 



expansion for {A[(l)i])^^fp.^y We also find from ([7]) that the higher-order terms in n give disorder correlation between 
j^.^j and the free energy W, which should be generally removed. We comment that substituting W for 
{A[<l)]) fj^^.^^ yields the moments of the free energy from the higher-order terms [sl. [l3l|. 



III. RENORMALIZATION GROUP 



The perturbation becomes more powerful combined with the renormalization group (RG), which consists of a 
coarse graining and a rescaling |14;] . The coarse graining means integrating over high-momentum components. Let 
the momentum space with a cutoff A he K. = {k : < \k\ < A}. Introducing i > 1, we divide K. into low- and 
high-momentum spaces defined as /C< = {fc : < |fc| < L~^A} and A^> = {k : L~^A < \k\ < A} respectively. The 
Hamiltonian H[(j);v] is also decomposed into H"^[(j);v] and where H^[(j)]v] contains only the low- momentum 

component (f>{p) {p G /C<), and H^[(j);v] denotes the remainders. Integrating over {q e /C>) in Z, we have the 
correction term 5H[(j), v] of the Hamiltonian generated as 

J n d0(<z)e-^'[^'"l -e-^^t*"!. (10) 

The rescaling procedure is carried out introducing the renormalized field 4>'{k) = (j){L~^k) and the disorder v' . 
Here, a number 9 and a renormalized disorder v' is determined in such a way that a main part of the Hamiltonian 
remains the same form: 



- H< [0; v] - 5H[4>; v] = -iJ [</>'; v']-5vq + -- - , (11) 

where 5vq means (/)-independent but w-dependent term, which should be added to the Hamiltonian because we put 
iJ[0;w] = 0. The dots in the right-hand side represents irrelevant terms other than 5vq, which can be neglected in 
low-energy physics. The same procedure leads to the RG transformation (RGT) for {A[(f)\) ^j^^.^^ [ll,[l3], where 
A[(f>\ contains low- momentum components of (j). Since the expectation value is independent of 5vq, the RGT becomes 

(^M).^[,;.]=(^[iV])^[,,„,]. (12) 

Taking the average over quenched disorder, we get 

E(^M)^[^;.]=E(A[iV'])^[^,„,]. (13) 



4 



Next, we calculate the counterpart in the replicated system. We apply to the RGT to Z in ([3]) before taking the 
random average: 



' 122 = 1 f^l<l>'a, ■■v']-nSvo 



which indicates that the renormalized replica Hamiltonian satisfies 

- H[cf)'; u'] = log E m't'>']~nSvo 
up to the irrelevant terms. Thus, the RGT for the replicated system is written as 



(14) 



(15) 



(16) 



Since there are no free replica sums in the right-hand side of (jlSp. u' explicitly depends on n through n5vQ. 
Furthermore, employing the cumulant expansion the renormalized coupling constants are expressed as polynomials 
in n. Combining ([9]), (fT6| and (fT3)) . we get 



(A[iV;]) 



n[(f>';u'' 



(17) 



It clearly shows that the right-hand side of (fT5|) can be computed from the RGT in the replicated system with the 
constant-term extraction. Note that Vq removes n dependence of u' . It implies that the beta function, which is 
defined by the linear response of u' under the infinitesimal change L ^ L + 6L, is independent of n. 



IV. EXAMPLE 



Now we exhibit a couple of examples. One of the simplest example is the random-field Gaussian model in d 
dimensions given by the following Hamiltonian 



U2(fcl,fe)'/>(fcl)0(fc2), 



where vi and V2 are quenched disorder. For simplicity, we ignore fluctuation of V2 and fix 

V2 (fcl,fc2) = {kf+t)fiK), 

where f{K) = {2Tr)'^S (K) with K = ^i, while vi obeys 

K(vi(fci),«i(A:2)) = A/(X). 



(18) 



(19) 



(20) 



The subscript k for the integral in means the measure d'^k/{2nY on K,. We regard the second term in p8)) as the 
unperturbed Hamiltonian and the first term as perturbation V . Since cumulants of vi higher than (j20p vanish, 
the higher order terms with m > 3 in the first line of ^ vanish when we take the disorder average. Hence we can 
readily derive the exact result 



E{c^{k^)c^{k2))H,{^.^,^ = (Go(fci) + AGo(fci)2)/(X), 



(21) 



where Go{k) = l/(fc^ + t). The same quantity is computed by the replica method as shown below. Employing ([5]) 
and (|20p . the replica Hamiltonian to (|18p is 



1 r 

Perturbative expansion with respect to A shows that 

(0l(fci)0i(fc2)}^^[0.„] 



(22) 



2 m 



Go(fci) + ^n^-^A^G^o+i(fci) I f{K). 



(23) 



FIG. 1: Vertices appeared in the Hamiltonian (|27|l . which are generated by k {V4) and k («3,ui) respectively. The open circle 
denotes vi. 




(a) (b) (c) 

FIG. 2: The one- loop correction discussed in the main text. A low momentum component is shown in a broken line, while the 
solid line means Go{q), q € /C>. A cross on a solid line depicts A carrying AG()(g)^, which is generated by k (vi,vi). 



We see that the constant term in n is identical with (|21|). as expected. On the other hand, if we take m ^ 00 in (j23|) . 
we obtain the exact two-point function in the replicated system: 

(0l(^l)'/'l(fc2))^^[0.„] 

= Go(fci) (l + UTTT^—ir^ f{K). (24) 



{kl+t- nA) 

It is argued in Ref.fToj that, when (j)'^ coupling constants are taken into account, the limiting procedure putting 
t = n A and then n — > in (|24p may generate a correction singular in n to the couphng constants due to n A in the 
denominator. However, it is easily checked that the higher-order terms with j > 1 in correspond to disorder 
correlation between {(t'{f^i)4'{k2)) H[(p-v] ^^'^ ^^^^ energy W, which should be removed. 

Next, in order to include interactions discussed above, consider the following Hamiltonian inspired by [Tol. fl^. : 

HW^^^^Y.n / vi{k^,...,ki)(ly{k^)---cj,{ki). (25) 

Z=l Jki,...,ki 

Here, the quenched disorder has the following non- vanishing cumulants in addition to (jl9p and (j20p : 

K(w4(fci,fc2,A:3,/c4)) ^Uif{K) 

K («3(fcl, k2,k3),V,{ki)) - -U2 f{K). (26) 

They yield the fohowing rephca Hamiltonian employing (|25p and (fS)). 

n 

u] = Hi [0; u]+ J2 (if'^"/' + ^t) ' (27) 

Q,/3 = l 

where Hi is given in ([^ and 



(7„;3= / /(if)0a(fcl)'^a(fc2)0a(fc3)0/3(A:4). (28) 

Here we focus on the one-loop corrections to ui having one A. We treat 7ii[</);w] as the unperturbed Hamiltonian, 
while the remaining terms in (j27p the perturbation graphically represented as Fig. [1] As we discussed in the above 
example, the free propagator is (j24p with nA removed. The n-dependent, one-loop diagrams having one A for u'l are 
presented in Fig[^ which are calculated as 



- (Suf + 6nuiU2 + Sn^Uj) / 

• Jq£K> 



Goiq)'. (29) 
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Vi Vi 





(b') 



FIG. 3: Diagrams corresponding to (b) and (c) in Fig. [5] before taking random average 



The novel limiting procedure proposed in 
/3i for gi = Aui: 



Pi 



[To'l is n ^ with g2 = nu2 fixed, yields the following beta function 

(6 - d)5i - (35? + 6Agi52 + i^^D- (30) 
In Ref. [lo| . a similar mechanism is argued in the context of failure of 



It means that 52 can affect flow of gi 
dimensional reduction. 

Contrary to this procedure, our general argument indicates that the second and the third terms in should 
be removed because they are generated by the disordered average in (|15p containing nSvQ. In fact, the diagrams in 
FiglJl which corresponds to those in Fig|5] before taking the random average, have nvi{q)Go{q)vi{—q) contained in 
perturbative expansion in uSvq. Thus the beta function for gi according to our argument becomes 



^1 = (6 - d)g, - 3gl 



(31) 



Here, 52 does not affect flow of gi within the one-loop correction, which leads to dimensional reduction. Although 
we have demonstrated the simpler model, a tedious but straightforward computation in the original model [10] shows 
that all the terms depending on n and remaining in the novel limiting procedure are caused by disorder correlations 
with nSvQ. These terms drastically change stability of fixed points and lead to breakdown of dimensional reduction 
[lo| . However, in the case when the higher order terms in n arc dropped in the beta function, this scenario docs not 
happen and dimensional reduction does occur il9J Although it is still debated whether dimensional reduction 
occurs or not near the upper critical dimensionsfsl. lid, fill [l2l . Iisl . ITol . I20I [2ll . [2^ . [2^ , the seeming ambiguity of the 
limiting procedure n ^ in perturbation cannot explain its breakdown. 



V. SUMMARY AND DISCUSSION 



We have shown that a physical quantity in a quenched disordered system calculated in a finite-order perturbation 
can be exactly derived by the replica method where the limit n ^ means to extract the constant term in n. 
The n-dependent terms give disorder correlations between the physical quantity we want to compute and the free 
energy, which should be removed. In this sense, the limiting procedure is uniquely determined as long as the physical 
quantity is represented in a polynomial in n. Our claim can provide a part of mathematical basis for various results 
by perturbative RG with the replica method where the limit is properly understood. 
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